Valency configuration of transition metal impurities in ZnO 
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Abstract 

We use the self-interaction corrected local spin-density approximation to investigate the ground 
state valency configuration of transition metal (TM=Mn, Co) impurities in n- and p-type ZnO. 
We find that in pure Zni_^TM^O, the localized TM^"'" configuration is energetically favoured over 
the itinerant d-electron configuration of the LSD picture. Our calculations indicate furthermore, 
that the (+/0) donor level is situated in the ZnO gap. Consequently, for n-type conditions, with 
the Fermi energy ep close to the conduction band minimum, TM remains in the 2+ charge state, 
whilst for p-type conditions, with ep close to the valence band maximum, the 3+ charge state is 
energetically preferred. In the latter scenario, modelled here by co-doping with N, the additional 
delocalized d-electron charge transfers into the hole states at the top of the valence band, and hole 
carriers will only exist, if the N concentration exceeds the TM impurity concentration. 

PACS numbers: 
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Based on the possible interplay between electronic properties and spin-functionality, di- 
luted magnetic semiconductors (DMS) are expected to play a major part in the development 
of the next generation of electronic devices.— The considerable challenge consists in designing 
materials that remain ferromagnetic above room temperature. In Mn doped GaAs, where 
ferromagnetism is well established, the currently highest achieved Curie temperature is 
~ 170 K. Ferromagnetism at even higher temperatures has been reported in some semicon- 
ductors (for a review, see for example Pearton et al."^), but doubts persist as to the carrier 
induced nature of the observed magnetic order. Recently, the prediction by Dietl et al^ 
of room temperature ferromagnetism in p-type Zni_a.Mn2.0 has generated considerable re- 
search activity, both in theory and experiment. So far, no conclusive experimental evidence 
has emerged that could either confirm or disprove the prediction. Various experimental 
investigations of the magnetic order in Zni.^^TMajO give contradictory results, ranging from 
spin glass^ and antiferromagnetic behaviour^'^ to ferromagnetism below- or even above^^^iiS 
room temperature. It has also been suggested that the observed ferromagnetism might be 
due to the formation of secondary manganese oxide phases.— Generally, the observed mag- 
netic properties of the DMS seem to depend strongly on the method by which they are 
synthesized. From all these conflicting results, the overall picture (though not universally 
accepted) seems to emerge, that Zni-j^TM^,.© will not be ferromagnetic,— unless additional 
hole or electron carriers are co-doped. Comprehensive overviews on the current state of 



research into ZnO spintronics can be found in references |l3lll4ljl5lll6 |. 

From the theory point of view, there are currently two main approaches to modelling 
dilute magnetic semiconductors that differ rather fundamentally in the way the transition 
metal (i-electrons are described. In the Zener model approach proposed by Dietl et al^ 
it is assumed that the on-site d-d correlations are significant enough to prevent hopping, 
causing the transition metal ^-electrons to remain localized on each their sites. In the band 
model approach it is assumed that the gain in kinetic energy associated with (i-electron 
delocalization is large enough to overcome the on-site correlations, and that consequently 
the (i-electrons fill up the resulting exchange split (i-band, situated in the case of ZnO in the 
middle of the gap. Both the Zener model description^^ and the band-electron pictur o^^i^^ 
have been put forward to predict high temperature (even room temperature) ferromagnetism 
in TM doped p-type ZnO. 

Despite this apparent agreement on ferromagnetism, the two models nevertheless differ 
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with respect to the proposed mechanism, which ultimately can be traced back to the very 
different assumptions regarding the ground state valency configuration of the TM impu- 
rities. With the self-interaction corrected (SIC)-LSD approximation, used in the present 
work,— i2S delocalized and localized d-electrons are treated on an equal footing, by removing 
an unphysical self-interaction, inherent in LSD, for every localized d-electron .— Thus the 
localized mean field scenario, and the delocalized bandstructure scenario can be compared, 
and the ground state configuration of the TM-ion can be deduced from the global energy 
minimum. 

I. SIC-LSD METHOD 

Whereas density functional theory (DFT) is in principle exact, the very popular local 
spin density approximation (LSD) used to describe exchange and correlation introduces an 
unphysical self-interaction. This can be demonstrated on a single electron moving in an 
external potential, which should be described by the total energy functional 

E""^^ K] = Ek^n [na\ + Kxi [na\ (1) 
but which in the LSD approximation is given by 

^^^^'^ [na\ = Ek^n [Uo] + Vext [flc] + Uh M + E^^"" [n„] . (2) 

Here, U^K] and E^f^[ Uq,] are respectively the self-Hartree and self-exchange-correlation 
energies, for the single electron density [uq], and together defining the self-interaction as 

€''' = UH[na\ + E^J''[n^]. (3) 

In the solid, using the SIC-LSD method, the self-interaction for each occupied orbital is 
subtracted from the LSD functional 

occ 

E'''^[m] = E^^^[n]-Y,5i^^^ (4) 

where the total energy E^^'~^[{^l)a}] is expressed as functional of the set of orbitals Tp^, and n is 
the total charge density of the system. For states of Bloch symmetry in an infinite crystal, i.e. 
band states, the self-interaction vanishes. To benefit from the self-interaction correction, a 
localized state must be formed. The ensuing many-electron wavefunction becomes a mixture 



of Bloch-type and Heitler-London type one particle states, i.e. delocalized and localized d- 
electrons coexist. The delocalized ci-states together with the remaining valence electrons 
continue to be treated in the LSD approximation, whilst those rf-states that are treated as 
localized gain the self-interaction (localization) energy, but loose the corresponding band 
formation energy. TM configurations with varying numbers of localized (i-states - in general 
different 'localization scenarios' - can be explored, and compared with respect to their total 
energies. Consequently, the preferred ground state configuration, as far as the number of 
localized d-electrons is concerned, can be determined from the global energy minimum. 

Different configurations of localized and delocalized d-electrons will give different nominal 
valencies, which in the SIC-LSD are defined as: N^ai = Z — Ncore — Nsic, where Z is the 
atomic charge, N^ore is the number of core and semicore electrons, and Nsic is the number 
of self- interact ion corrected (i-electrons respectively. According to this definition, a divalent 
TM valency refers to the localized configuration of the Mn ion (Mn^+) , and the localized 
(f Co ion (Co^"*") respectively. The trivalent configuration has one less (i-electron localized, 
and thus refers to the localized (i^ Mn ion (Mn^"*") and Co ion (Co^+). The LSD picture 
corresponds to the scenarios with all (i-electrons delocalized, i.e. the d^ configuration for the 
Mn ion (Mn^+) and the Co ion (Co^+). 

II. BAND PICTURE VERSUS LOCALIZED PICTURE 

We use the SIC-LSD method to study the validity of the two above mentioned electronic 
structure models of TM impurities in ZnO, and to determine which of the two assumes 
the correct TM ground state configuration, based on total energy considerations. When 
doping ZnO with Mn/Co, the TM ions occupy the Zn sites without changing the wurtzite 
structure. In our calculations, Zni^^jTM^O is modelled, by substituting a single Zn in 
a (2x2x2) supercell consisting of 16 ZnO formula units, with a TM impurity. The SIC- 
LSD scheme is implemented within the all-electron tight-binding linear-muffin-tin orbitals 
(TB-LMTO) method, with the atomic sphere approximation (ASA). Spin-orbit coupling is 
fully included in the self-consistency cycles. Empty spheres are introduced on high symmetry 
interstitial sites. The electron charge density was sampled using 27 A;-points in the irreducible 
wedge of the Brillouin zone. All calculations use the experimental lattice parameter. The 
resulting densities of states (DOS) for the delocalized d-electron scenario (LSD picture), and 
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FIG. 1: Total (black), and Mn-d projected (red) DOS as a function of energy, in states per eV, 
of Zni_i.Mn2,.0 a) in the delocalized (band picture) Mn^"*" configuration and b) in the localized 
divalent Mn'^^ configuration. 

the localized scenario are shown in Fig. and Fig. CJd respectively. 

In Fig. Ht., the delocalized electrons fully occupy the spin-up channel of the exchange 
split (i-band whilst the spin-down band remains empty. The proponents of the band pic- 
ture conclude that ferromagnetism, due to double-exchange, can not occur in this scenario. 
The same argument however is less useful for explaining the absence of ferromagnetism in 
Zni_j;Coa;0, where the two additional d-electrons start filling up the spin-down band, and 
the corresponding LSD DOS is reminiscent of a half-metal. In Fig. ^p, the 5 Mn d-electrons 
have been Sl-corrected, and are therefore no longer available for band formation. Instead 
they form localized levels at energies below the valence band maximum (VBM). Note that 
the SIC-LSD, like any other density-functional theory, does not provide excited state prop- 
erties. In particular, the approach cannot accurately determine the removal energies of 
localized states, and the bare Sl-corrected d-bands always appear at too high binding ener- 
gies due to the neglect of screening and relaxation effects^-^ and are therefore not shown in 
the DOS plots. The Fermi level is situated in the middle of the ZnO gap, and there are no 
free carriers, which for the proponents of the localized model explains the absence of carrier 
mediated ferromagnetism. The corresponding DOS for Zni_j.Co^O is quite similar to Fig. 



^P, but with seven Co (i-states now Sl-corrected in the locahzed scenario. Comparing the 
total energies for both models (columns 2 and 5 in Table I), we find that the localized Mn^+ 
and Co^^ scenarios are energetically more favourable than the corresponding band scenarios, 
by respectively 3 eV for Zni_^.Mn2.0 (row 2) and 5.5 eV for Zni_a;Coa.O (row 8). 

The band picture is clearly not preferred energetically, relative to the localized picture, 
however we still need to investigate whether the global energy minimum is obtained for an 
intermediate valency configuration of coexisting localized and delocalized rf-electrons. In 
Fig. 121 the DOS of Zni_2.Mna;0 is shown with the Mn ion in the Mn^+ configuration, i.e. 
only four of the d-electrons are localized by the SIC. The remaining d-electron, which is 
now treated in the LSD approximation, fills up the small d-peak situated in the middle 
of the ZnO gap. From Table I (row 2, columns 2 and 3), we see that this configuration 
is energetically less favourable than the TM^"*" configuration. Similarly for Co doped ZnO 
the global energy minimum is obtained in the Co^+ scenario (Table I, row 8), and we 
conclude that both Zni_a,Mna,0 and Zni_2^Coa;0 prefer the localized rf-electron ground state 
configuration. The TM^+ configuration of transition metal impurities in ZnO is in agreement 
with experiment,— and in support of the electronic structure underlying the Zener model, 
rather than the standard band picture. 



TABLE I: Total energies in eV of Mn/Co doped p- and n-type ZnO. 



Compound 


Mn2+(d5) 


Mn3+(d4) 


Mn4+(d3) 


Mn7+(d0) 


Znis/ieMni/igO 


-113.13 


-112.61 


-112.03 


-110.12 


Znis/igMni/ieOis/igNi/ie 


-111.40 


-113.70 


-112.58 




Zni5/i6Mni/i60l4/16N2/16 


-110.80 


-111.48 


-111.35 




Zn^/igMni/igGai/ieO 


-115.91 


-115.40 






Zni4/i6Mni/i6Sni/i60 


-110.49 


-109.89 








Co2+(d7) 


Co3+(d6) 




Co9+((i0) 


Ziiis/ieCoi/ieO 


-114.04 


-113.40 


-112.70 


-108.38 


Znis/ieCoi/igOis/ieNi/ie 


-112.31 


-113.90 


-113.25 




Zni5/i6Coi/i60i4/i6N2/i6 


-110.28 


-112.01 


-111.79 




ZnM/ieCoi/igSni/ieO 


-111.48 


-110.89 







fi 



10 







-10 





' 1 ' 1 ' 1 ' 


1 1 1 1 1 1 1 1 1 1 

i ZnO:Mn(d4) 

I 

i 

i 

I 




. 1 . 1 . 1 . 1 ' . 1 . 1 . 1 . 1 . 



-4 -3 -2 -1 1 2 3 4 5 6 7 8 

FIG. 2: Total (black), and Mn-d projected (red) DOS, in states per eV, as a function of energy, in 
eV, of Zni_2.Mna;0 in the trivalent Mn^+ configuration. 

In Fig. 121 the gap state that forms when treating an additional (i-electron as band electron, 
is relatively narrow, and completely filled. The corresponding gain in band formation energy 
is insufficient to overcome the loss in Sl-energy, and the (i-electron prefers to remain localized. 
If however the Fermi level could be positioned below this gap state, for example through 
additional acceptor dopants, one could very well imagine that the resulting gain in energy 
from charge transfer to the Fermi level would be large enough for the delocalization to occur. 
The dependence on Fermi level position can be expressed through the formation energy 

EfiZnO : ^M^ q) = Et^tiZnO : TM\ q) - Etot{ZnO) - f^zn + I^tm + qep (5) 

where (ZnO:TM'',q) refers to an impurity of valency v, and charge q, fizn and fiTM are the 
respective chemical potentials, ep is the Fermi energy, and Etot the total energy of the system 
with impurity. Three scenarios need to be considered: In the neutral charge states, the sub- 
stitutional impurity can either assume the divalent configuration, TM^"*" = (ZnO:TM^"'",0), 
which does not introduce states in the gap (Fig. Q), or the trivalent configuration TM^"*" = 
(ZnO:TM^"'",0) which has an electron rf-state in the gap (Fig. |2I). In the charged scenario 
(ZnO:TM^"'",+) the delocalized (/-electron has been transferred to the conduction or acceptor 
band, ionizing the TM impurity into the positive charge state. The formation energy as a 
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ZnO:TM 




Fermi energy ^BM 

FIG. 3: Schematic plot of the formation energy as a function of Fermi energy for ZnO:TM, illus- 
trating the link between the localization and delocalization picture and the TM donor ionization 
level in the band-gap. Only the (+/0) level between the positively charged TM^^ and the neutral 
TM^+ valence configurations is indicated. 

function of Fermi level is plotted in Fig. |21 The two horizontal lines correspond to the 
neutral scenarios, that are independent of the Fermi level position. The formation energy 
for the charged scenario changes linearly with Fermi level, as follows from equation (5), and 
as indicated by the skew line. The delocalization process 

rM^+ TM^+ + electron,^ (6) 

will go ahead, if the position of the Fermi energy is such that Ej(ZnO:TM^"'",+) < 
Ej(ZnO:TM^"'",0). The donor transition level (+/0) is defined by the Fermi energy above 
which the TM impurity is in the divalent state, and below which it is in the positively 
charged trivalent state. From SIC-LSD calculations, we have determined that the donor 
level in Zni_2:TMa;0 is positioned more or less centrally in the ZnO gap^. The same con- 
clusion can be drawn qualitatively from applying the numbers of Table I to Fig. El Since 
the neutral trivalent state (ZnO:TM^"'",0) has a gap state in the high energy part of the 
band gap, and since it is situated only ~ 0.5 eV above the neutral divalent (ZnO:TM^+,0) 
configuration, the donor level (+/0) will be situated in the 3.4 eV wide ZnO gap, as is also 
found experimentally.-^ 

To investigate the actual effect of additional dopants on the electronic structure of 
Zui^^TMa-O, we use the SIC-LSD to determine the ground state configuration for both 
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FIG. 4: Total (black), N-p(blue), and Mn-d(red) DOS as a function of energy, in states per eV, of 
Zn;^5/X6Mn;^/;^gOi5/igN;^/;^g [ii tLc trlvalcnt Mn^+ ground state configuration. 

p- and n-type scenarios. Acceptor (p-type) doping is achieved by replacing a single O in 
the supercell by N, whereas n-type doping is modelled by replacing a single Zn atom by 
either Ga or Sn. The total energy results for the different dopants and different TM valency 
configurations are given in Table I. We find that, with N as additional dopant, the trivalent 
TM configuration becomes energetically most favourable, both for Mn (row 3) and Co (row 
9). Nitrogen, having one p-electron less than O, acts as acceptor when introduced into ZnO, 
with the corresponding acceptor level situated near the top of the valence band (Fig. |^. 
Given that the (+/0) level is situated above the acceptor level, the delocalization process 
(6) can go ahead, leaving behind an empty d-band (gap state). The energy gain from charge 
transfer and hybridization in the TM^"*" configuration is obviously large enough to overcome 
the corresponding loss in SIC (localization) energy. 

In their paper on the possibility of magnetism in Mn doped ZnO, Dietl et al."^- state 
that with n-type doping one can not expect ferromagnetic order, except maybe for very 
low temperature, but they predict room temperature ferromagnetism in the p-type scenario, 
where the hole carriers introduced by co-doping with N will mediate the interaction. The 
assumption is that the ground state configuration remains Mn^"*", based on the argumen- 
tation that the (+/0) level is situated below the VBM.^^ This picture is thus in contrast 
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to our calculation, which finds the (+/0) level situated above the N based acceptor states, 
and results in the Mn^"*" ground state configuration. The importance of determining the 
correct ground state configuration becomes clear from Fig. where the actual trivalent 
configuration implies charge transfer into the acceptor states, and absence of hole carriers, 
while a hypothetical divalent scenario with five localized d-states would imply absence of 
the gap state and charge transfer, leaving the N acceptor states partially unfilled (Dietl's 
scenario). From the SIC-LSD calculations it follows that, the band description of the d- 
states does not fully account for the correct electronic structure of the Zni_j.TM^O ground 
state, even when co-doped with N. But also the Zener model, with hole carriers mediating 
the magnetic interaction between localized spins residing on TM^+ ions, is not a true repre- 
sentation of the ground state. There are no hole carriers in either Znis/igMni/ieOis/ieNi/ig 
or Znis/igCoi/ieOis/ieNi/ie and carrier mediated ferromagnetism can therefore not occur. 

It seems plausible that increasing the concentrations of the dopants, [TM] and [N], will 
not change the overall picture of fully compensated acceptor states, as long as [TM] = [N]. 
However, the qualitative picture changes considerably if we increase the relative amount of 
N impurities, i.e., if [N]>[TM]. Substituting two of the O atoms by N, but with a single TM 
impurity in the 32 atom ZnO supercell, we find that, contrary to what one might expect, 
overdoping with N does not result in a further delocalization, i.e., a transition from TM^"*" 
— > TM^+ + e~. From Table I we see that the TM^+ configuration gives the lowest total 
energy, both for Zni5/i6Mni/i60i4/i6N2/i6 (row 4), and Zni5/i6Coi/i60i4/i6N2/i6 (row 10). 
The additional N impurity is thus not compensated by TM c/-electrons which results in an 
impurity band at the top of the valence band, which is now only partially filled. 

Co-doping with N has revealed itself to be rather difficult,— which according to calcula- 
tions by Yamamoto et alM can be related to the corresponding increase in the Madelung 
energy. Co-doping with electron donors, on the other hand, can be readily achieved. The 
theoretical prediction, that high temperature ferromagnetism will not occur in n-type Mn- 
doped ZnO, seems to be confirmed by experiment with respectively Al and Ga as co-dopants, 
however ferromagnetism has been reported in Zni_a;Mni,0, when co-doped with Sn.^^ In the 
case of n-type doping, the substitution of Zn by either Ga or Sn introduces respectively one 
or two additional electrons into ZnO, which results in an impurity band situated below the 
conduction band minimum (Fig. Ej). From Table I, rows 5, 6, and 11, we see that in this 
case the divalent TM configuration remains energetically more favourable, i.e. given that 
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FIG. 5: Total DOS, in states per eV, as a function of energy, in eV, of a) Zni4/igMni/igGai/i60 
and b) Znx4/igMni/i6Sni/i60, in the divalent Mn^^ configuration. The green, red, and black lines 
respectively represent the 0-p projected, Mn-d projected, and total DOS's. The Fermi level (zero 
energy) is now situated in the conduction bands. 

the Fermi level is here situated above the (+/0) level, delocalization does not occur, and 
respectively five (Mn) or seven (Co) d-electrons remain localized at each of the respective 
sites. With or without n-type co-dopant the divalent state (of both Mn and Co) invariably 
is 0.5-0.6 eV more favorable than the trivalent one, i.e. the n-type doping does not infiuence 
the TM configuration at all. 

With respect to the calculated magnetic moments of the TM impurities in the vari- 
ous scenarios and valency configurations, we find that the orbital moment remains largely 
quenched. With no additional co-dopants, i.e. rows 2 and 7 of table I, the spin moments 
are approximately 4.5 /i^ and 2.7 /i^ for Mn and Co impurities respectively. These mo- 
ments remain virtually unchanged regardless of the considered configuration, indicating for 
example that the SIC localized (i-state of the Mn^"*" configuration has the same moment 
as the filled narrow d-band state of the corresponding Mn^"*" configuration. This situation 
changes when N is codoped (rows 3 and 9 of Table I). Here in in the TM^"*" scenario, the 
moments are respectively 4.35 fiB for Mn in Znis/igMni/ieOi/igNi/ie and 2.53 /i^ for Co in 
Znis/ieCoi/ieOi/igNi/ie. However for the corresponding TM'^^ configuration, a (i-electron 
has charge transfered into the hole state, with the moment being determined by the remain- 
ing SIC localized d-states, i.e. Mn((i^) and Co{d^) . As a result, the magnetic moment of 
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Mn decreases to 3.7 /^s, and the magnetic moment of Co increases to 3.0 ij,b- 

III. CONCLUSION 

In summary, we have studied the electronic structure and different valency configurations 
of Co and Mn impurities in p- and n-type ZnO using the SIC-LSD ab initio method. Prom 
total energy considerations we find that the TM d-states remain locahzed if no additional 
hole donors are present. The TM^+ remains the ground state configuration under n-type 
conditions. The TM^+ becomes more favourable in p-type ZnO, and the holes required for 
carrier mediated ferromagnetism occur only when the N co-dopant concentration exceeds 
the concentration of the TM impurities. 
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